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ABSTRACT 

We have derived a new expression for the thermohaline mixing coefficient in 
stars, including the effects of radiative levitation and external turbulence, by 
solving Boussinesq equations in a quasi-incompressible fluid with a linear ap- 
proximation. It is well known that radiative levitation of individual elements 
can lead to their accumulation in specific stellar layers. In some cases, it can 
induce important effects on the stellar structure. Here we confirm that this ac- 
cumulation is moderated by thermohaline convection due to the resulting inverse 
/i-gradient. The new coefficient that we have derived shows that the effect of 
radiative accelerations on the thermohaline instability itself is small. This effect 
must however be checked in all computations. We also confirm that the pres- 
ence of large horizontal turbulence can reduce or even suppress the thermohaline 
convection. These results are important as they concern all the cases of heavy el- 
ement accumulation in stars. The computations of radiative diffusion have to be 
revisited including thermohaline convection and its consequences. It may be one 
of the basic reasons for the fact that the observed abundances are always smaller 
than those predicted by pure atomic diffusion. In any case, these processes have 
to compete with rot at ion- induced mixing, but this competition is more complex 
than previously thought due to their mutual interaction. 

Subject headings: stars: abundances, stars: hydrodynamics, stars: convection, 
stars: diffusion, Chemically peculiar stars 
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1. Introduction 

During many decades, the theory of thermohaline convection was developed in the frame- 
work of oceanography, as it represents on e of the majo r physical factor s leading to the large 
stream s of the earth ocean (jSternl ( J196CH ) . ISternl (Il967l ). I Turner! ( 119731 ). iGargett k, Ruddick 
( 2003h ). It was not much applied t o astrophysics, except for very few papers, like Ulrich 
and iKippenhahn et al.l (119801 ) (hereafter KRT80). At the present time, thermohaline 
convection is recognized as a major mixing process which occurs in stars in the presence 
of an unstable //-gradient associated with a stable temperature gradient, provided that the 
//-gradient is not large enough to drive dynamical convection. If a blob of stellar material be- 
gins to move down, it exchanges heat with its surr oundings more rapidly than particles, and 
goes on falling, thereby triggering the instability ( IVauclairl (120041 ) and references therein). 



Various situations lead to thermohaline convection in stellar interiors. They may be 
classified as follows: 



Accretion of heavy matter onto the star, which may come from planetary material 
( Vauclair 2004 ; Garaud 2011 ; Theado Sz Vauclair 2012h or from an evolved companion 
in a binary system ( Stothers fc Simon 19691 ; Stancliffe et al. 2007 ; Thompson et al. 
2008h . 



Local reductio n of the mean molecular weight due to nuclear reactions, for example 



in Red Giants (Eggleton. Dearborn fc Lattanzio 2006, 



2008 



Charbonnel k Zahn 2007 



Charbonnel fc Lagarde 



2010 



Stancliffe 



2010 



Denissenkov 



201 



Heavy elements accumulation induced by atomic diffusion (ITheado et al.l 120091 ) or he- 
lium accu mulation induce d by diffusion in a mass loss flux like in main-sequence helium 
rich stars ( IVauclairl 1 19751 ) . 



The variations with depth of the radiative accelerations on specific elements can lead 
to their accumulation or depletion in various layers inside the stars, which may have stron g 
consequences on the stellar structure and evolution (IRicher et al.ll2000b iRichard et al.ll200ll ). 
The accumulations of iron and nickel, which represent important contributors to the opac- 
ity in some stellar layers, may lead to extra convective zones, and ca n even, in some 
cases, trigger stellar pulsations through the iron-ind uced K-mechanism (jCharpinet et al. 



19961 ; iPamyatnykh et al.ll2004j ; iBourge fc Alecianll2006l ). However, in all these previous stud- 
ies, the thermohaline convection induced by the heavy elem ents accumula t ion w as not taken 
into account. First computations of this effect were given in lTheado et al.l (120091 ). who found 
that the accumulation of heavy elements was attenuated by the thermohaline mixing pro- 
cess, but not completely suppressed. The helium settling which occurs simultaneously with 
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the heavy element accumulation induces a stabilizing contribution to the global /i-gradient, 
which helps keeping some accumulation. 



Recent n umerical simulations of thermohaline convection flDenissenkovl (120101 ). here- 
after D10, and lTraxler et al.l (120111). hereafter TGS11 ) brou ght a new light on these previous 
computations. As discussed by iTheado fc Vauclair ( 2012 ). the thermohaline mixing coef- 
ficients previously u sed in the litterature have to be revised. Also, in the computations by 
Theado et al.l ( 120091 ). the effect of the radiative acceleration on the falling blobs was not 
introduced. Finally t he influence on the thermohal i ne con vection of an external turbulence, 
already addressed by iDenissenkov fc Pinsonneaultl (120081 ). is derived and introduced in the 
present context. 

In the present paper, we analyse this situation and we give a new derivation of the 
thermohaline mixing coefficient, including radiative accelerations and/or external turbulence, 
consistent with the numerical simulations. This new coefficient will be used in the future to 
compute the evolution of the element abundances inside stars, including the simultaneous 
effects of atomic diffusion and thermohaline mixing. The introduction of thermohaline mixing 
in the computations may have important consequences on the stellar structure and evolution. 



2. Modelling thermohaline convection inside stars 

Thermohaline convection is not a real diffusion process. It is a special kind of convec- 
tion which presents, at least in the ocean case, elongated cells called "salt fingers". This 
convection, induced by an unstable chemical stratification, leads to chemical mixing. For 
this reason, and for simplicity, it is generally modelled in the same way as diffusion, with 
a mixing coefficient whose expression has to be deduced from physical considerations. In 
the past, the mixing coefficients proposed for stellar conditions were ve ry unc e rtain, with 



unknown normalisation values. For example, the expressions given by lUlrichl (ll972f ) and 
KRT80 differ by nearly two orders of magnitude. Recently, numerical simulations like those 
of TGS11 have led to real improvements, as the mixing coefficient must account for the 
chemical fluxes obtained from the simulations. 

An important physical parameter in the study of thermohaline convection is the so- 
called "density ratio" in oceanography, which may be transposed to the astrophysical case 
as: 

R« = — T^-T— (1) 

where V a d and V are the usual adiabatic and local (radiative) gradients d In T/d In P, and 
= d In fi/d In P. 
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As discussed in I Vauclairl (120041 ). in the absence of external turbulence or radiative levi- 
tation, thermohaline convection may develop only if: 



1 < R n < 



(2) 



where r is the inverse Lewis number, equal to the ratio of the particle diffusivity to the 
thermal diffusivity. 

These two boundaries, which are well known in oceanography, have important physical 
meanings. For thermohaline convection to occur, the medium must be stable according 
to the Ledoux criterion (no dynamical convection), but the relative variations of T and 
\i inside a falling blob must be large enough, compared to the surroundings, to trigger the 
instability. This last condition imposes that the ratio of the thermal diffusivity to the particle 
diffusivity, which induces the variations of these two quantities inside the blob, be larger than 
the external ratio of the temperature gradient to the particle gradient. 

If Rq is smaller than one, the /i-gradient triggers dynamical convection. If Rq is larger 
than 1/r, the initial /i-gradient is not large enough to overcome the stabilizing effect of the 
temperature gradient. This means that the expression of the thermohaline coefficient must 
behave in the right way at these two limits: it must become very large for R close to one, and 
it must vanish for Rq close to 1/r. This behavior is not reproduced by any of the previously 
used thermohaline mixing coefficients, except for that proposed by D10, as discussed below. 



Ulrichl ( 119721 ) and KRT80 proposed to model the thermohaline mixing process in terms 



of a diffusion coefficient propor tional to the inverse /i-gradient, which can simply be written 



see 



Theado fc Vauclair I (l2012h ) 



D th = C t k T R^ 



(3) 



where kx is the thermal diffusivity, and Ct an unknown factor including the aspect r atio 
(length/width) of the fingers. C t is ev aluated as 12 by KRT80, an d 658 by lUlrichl ( 119721 ). In 
their computations of red giant stars, ICharbonnel &: Zahnl (120071 ) considered this parameter 
as free and adjusted it as C t = 1000. Note that this expression does not satisfy any of the 
two physical limits: D t h goes to infinity for Rq = and not Ro = 1 as needed, and it does 
not vanish for R = 1/r. 

The recent numerical simulations by TGS11 considerably improved our knowledge of 
thermohaline convection. They found universal scaling laws for the turbulent heat and 
compositional fluxes, valid for different sets of Prandt and Lewis numbers. Fitting their 
numerical results by an empirical law, TGS11 gave an expression for D t h which vanishes 
for Rq = 1/r, but reaches a finite limit for Rq = 1, which is not correct. Interestingly 
enough, a close look at their simulation points (their Figure 2) shows that the heat and 



-5 - 



compositional fluxes do increase for small Ro, far above their empirical relation, suggesting 
that this relation should be revised. 



Denissenkovl ( 120101 ) (D10) derived a thermohaline mixing coefficient using a straight- 



forward linear treatment of Boussinesq equations in a quasi-co mpressi b le flu id. He found 



an expression which nicely satisfies the boundary conditions. As lUlrichl ( 119721 ) and KRT80, 
he introduced an adjustable parameter, C t , depending on the shape of the convective cells 
(aspect ratio of the fingers) and tried to constrain it by 2D simulations. Here we prefer using 
the more recent 3D simulations of TGS11 to adjust this parameter. We find that the D10 
expression nicely fits the simulations, including the increase of D t h when Rq becomes close 
to one, which is visible on the simulation but was not introduced in the TGS11 expression: 

D th = C t k T (Ro - 1) -1 (1 - R t) (4) 

Figure 1 presents the behavior of all these coefficients as a function of Rq, in stellar 
conditions. As an example, they have been computed for a 1.7M Q , 403 Myr model. Here the 
1/r value, which is related to the rapid decrease of D t h for large values of -Rp, is 6.8.10 8 , which 



i s larg er than the values generally given in the literature (e.g. TGS11 or lTheado fc Vauclair 



( 120121 ) ) . The basic reason for this difference is that this parameter is generally computed 
with helium as the diffusing element, whereas here we take into account the fact that the 
diffusive element is iron, which diffuses much more slowly than helium (C.f. Figure 2). From 
all the considerations discussed above, the coefficient proposed by D10 is clearly the best 
one, as it satisfies the results of the 3D simulations, as well as the boundary conditions. 
However, this coefficient does not take into account external turbulence, neither the effect of 
the radiative acceleration on the falling blobs, which may be important when the instability 
is induced by heavy element accumulation due to radiative levitation. In the following, we 
derive a new thermohaline coefficient in the same way as D10, including these two effects. 



3. New derivation of the thermohaline mixing coefficient including radiative 

acceleration and/or external turbulence 

In this section we give a linear analysis of thermohaline convection in the presence 
of atomic diffusion and/or external turbulence. We begin with Boussinesq equations that 
describe motions in a nearly incompressible stratified viscous fluid, in the same way as 
D10, except that we add in the equations the terms due to atomic diffusion and turbulent 
diffusivity. 

Atomic diffusion is a selective process related to the fact that the stellar gas is composed 
of many different elements, each of them behaving its own way in the stellar gravity and 



2xl0 8 



4x10 s 



6xl0 8 



Fig. 1. — Examples of various prescriptions for the thermohaline diffusion coefficients as 
computed in a 1.7M , 403 Myr model, with a f/r value of 6.8.10 s . In the main frame, the 
dotted cu rve represents the KRT80 co efficient, with C t = 12, as a function of Ro, the dashed 
curve the ICharbonnel &: Zahnl ( 120071 ) coefficient, with C t = 1000, the dotted-dashed curve 
represents the TGS11 expression and the solid curve the D10 coefficient, adjuste here with 
C t = 120 to fit the TGS11 simulations. In the small frame, the solid curve is a zoom of the 
D10 coefficient for Ro < 6.10 4 , the two other curves are obtained for horizontal turbulent 
diffusion coefficients equal to 3000 cm 2 .s _1 (dashed curve) and 15000 cm 2 .s _1 (dotted-dashed 
curve) . 
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radiation fields. If the radiative acceleration gnii) acting on a given element (say iron) 
is locally larger than gravity, this element moves up with respect to the average stellar 
matter (mainly composed of protons). They move individually until they share their acquired 
momentum with the surroundings, which occurs after a collision time scale. The overall effect 
of this individual radiative acceleration on the medium is of order gn = T,iXign(i), where 
Xi is the mass fraction of the considered element. The global radiative acceleration on the 
medium g^ is generally several orders of magnitude smaller than the gravity, so that the 
medium as a whole is not destabilized (Figure 3). However, if the radiative acceleration on 
the considered element decreases with increasing radius, the diffusion flux decreases as well 
and this element accumulates in some layers inside the star. As a result, the internal stellar 
structure may be modified, due to changes in the opacity. Meanwhile, an inverse /i-gradient 
is created, which may become unstable against thermohaline convection. 

External turbulence (for example rotation -induced mixing) may a l so coe xist with ther- 



mohaline convection. In previous works (e.g. ICharbonnel fc Lagardd ( 120101 )). such a situ 



ation was treated by simply adding the m ixing coefficients of the two process es, computed 



separately, in spite of the analysis given by iDenissenkov fc Pinsonneault 



(120081 ). who showed 



that the thermohaline mixing coefficient must be modified in this case, due to the action of 
the external turbulence on the falling blobs. Here we recalculate and confirm this effect and 
introduce it in the final expression of the thermohaline mixing coefficient. 

Following D10, we first write the Boussinesq equations in a nearly incompressible fluid 
in the following form: 

^ + (v, V)v = ^g e + <5g e + z/V 2 v (5) 
8T 

— + (v, V)T = k* T V 2 T (6) 

^ + (v,V) / u = A;;[V 2 / i + /(ge,r)] (7) 

In the first equation, we have to take into account the fact that the gravity acting on 
the fluid is now an effective gravity defined as g e = g — g^, where g^ is the overall radiative 
acceleration acting on the medium, namely g^ = EjX^^(i). As already mentionned, g^ is 
most generally negligible compared to g. However, when the effective gravities inside and 
outside the blobs are substracted, the gn terms are the only remaining ones, as they are 
different in the falling matter and in the surroundings. 

The second equation, which describes the heat transfer between the falling blobs and the 
external medium, remains unchanged except that, in case of external horizontal turbulence, 
the heat transfer induced by the extra mixing must be taken into account in the thermal 
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diffusivity. When the Peclet number, ratio of the horizontal turbulent diffusion coefficient 
to the thermal diffusivity is much smaller than one, k^ reduces to kjn, which is given by: 

4acT 3 



2>KCpp 2 



But when the Peclet number becomes larger than one, the horizontal diffusion coefficient 
must be added to kx in the expression of k^. 

The third equation concerns the diffusion of particles between the blobs and their sur- 
roundings. The coefficient k^ is the so-called "particle diffusivity" which, in hydrodynamical 
contexts, does not take into account the fact that stars are made of multicomponent gases. 
It is generally computed in an approximate way, and the values cited in hydrodynamical 
papers often correspond to the diffusion coefficient of helium in a hydrogen rich medium. 
Sometimes, it is even identified to the molecular viscosity, which is nothing else than the 
hydrogen- hydrogen self diffusion coefficient. In the atomic diffusio n context, the d iffusio n 



coefficients for each element are computed in a consistent way (e.g. iPaquette et al.l (119861 )). 
It is not possible, in the present situation, to take into account the separate diffusion of all 
the elements through the blobs boundaries, but one must at least be careful to use for k^ the 
diffusion coefficients corresponding to the most important elements (e.g. iron). The ratio of 
the diffusion coefficients for two different ions (1) and (2) of respective charges Z\ and Z 2 is 
approximately (Z 2 /Zi) 2 . Using that of iron leads to a smaller diffusivity than using that of 
helium, which would not be correct in this case (Figure 2). 

In a general way, the particle diffusion between the interior and the exterior of the 
falling blobs is not only due to the abundance gradients, but also to the effects of gravity, 
radiative acceleration, tem perature gradient, etc., which occur in the vertical direction (e.g. 



Vauclair fc Vauclairl (J1982J)). This is added as /(g e ,T) in Equation 7. However, as pointed 
out by previous authors (e.g. D10), the vertical selective diffusion velocity is negligible 
compared to the downward velocity of the blobs so that we can ignore this effect. On the 
other hand, as mentionned above, the average effect of the radiative acceleration on the 
medium has to be taken into account in the computation of the velocity (Equation 5). 

Finally, in case of external turbulence, the mixing of particles between the blobs and 
the inter-blobs is increased. The horizontal component of the extra turbulent diffusivity has 
then to be added to fc M in equation 7. We have written the total diffusivity as A;*, which 
reduces to k^ when no external turbulence is involved. 

In the following, we use the cartesian coordinate system with the vertical z axis pointing 
in the direction opposite to the gravitational acceleration, the horizontal x axis in a star's 
meridional plane and the horizontal y axis perpendicular to them. Linearizing equations 
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(5, 6, 7), with w representing the vertical component of the blob velocity, we obtain the 
following system of equations, similar to equations 5 to 7 of D 10 except for the gravity term 
and the effect of external turbulence included in and A;*: 

= g e (a5T - 06 + uV 2 w + 5g e (9) 

OSli dfi 2 

where a and j3 are the usual coefficients of the equation of state: 

?£ = -adT + pdn (12) 
Po 

Note that for simplicity we keep here the notation dT/dz but we have to remember that in 
stars, due to the temperature stratification, this is equivalent to (V ra d — V ' a d)(T / Hp). 

When a heavy blob falls down, the conditions inside it vary in a non adiabatic way, 
due to heat and particle diffusion with the surroundings. The vertical variations of T and \x 
inside and outside the blob may be expressed in the following way: 

9T ^d-W (13) 



dz dz 
dji dpLo 



i-<y (14) 



dz dz 

where the coefficients 5t and 5^ have values between (equality between the internal and 
external values of T and /x) and 1 (adiabaticity, no internal variations of T and //). They 
may be defined from the following relations which give the differences ST and 5 fx between 
the inside and the outside of the blobs. 

^=-if«r (15) 

where T and /x are the temperature and mean molecular weight outside the blobs and t is 
the unknown length of the elongated cells (fingers). 

The blobs stop falling when the buoyancy inside and outside are identical, that is when 
pg e are the same inside and outside. 
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As in D10, we transform equations 9, 10, 11 in terms of w, St and 5 M and we search 
for solutions in the form of exp(cxt) x exp(i[k x x + k y y}). But before that we have to express 
the new term Sg e in the same way (in terms of St and S^ ). In the expression of the 
effective gravity, the only term which is different inside and outside the blobs is the radiative 
acceleration, which is a function of both the temperature and the chemical composition. 
Although these two variables are not completely disconnected, we will separate them in first 
approximation. This approximation will be discussed in the following section. We write: 



or 

fdg e dT £ dg e d/i i,\ fidfi dg e T dg e \ 

In the following, for simplicity we assume a perfect gas law, so that a — 1/T and 
(3 = l//i. Also, for simplicity we write here kT and k^ instead of k T and k* Defining g e ^ as 
<91ng e /<91nT and g e ^ as d\ng e /d\n/i, the preceding equations may then be reduced to : 

4cr 3 + 2(k T + k ll + 2u)k 2 a 2 + 

{[k T k, + 2v{k T + k„)]* 4 + 2gj^{R Q - 1)} a 

-g e — (k T - R k^k +h>k T k tl ¥ > 

din /j,o , 2 din /j,o , 

+2ag e — {g e>lt + Rog e ,T> + k g e d (k T g e ,n + Rok^g e ,T) = (19) 

with k 2 = k 2 x + k 2 y 

Taking into account that the thermal diffusivity is much larger than the particle diffu- 
sivity and the molecular viscosity, we neglect the same terms as D10 and this equation is 
simplified in: 



= k2 ( k T ~ Rpkf,) + (k T g e ^ + Rpk^g e)T ) , . 

°~ 2 (R -l)-(g e4i + R g eT ) [ > 

The thermohaline mixing coefficient, of the order of l/2a£ 2 , becomes: 

D = Ck Ai-Ror*n(g e ,» + R T*g t) 

T (Ro -1)- tie* + Roge,T) 

where kT has been replaced by k T and the inverse Lewis number r has been replaced by 
r* = k*Jk T . C t is the usual constant related to the aspect ratio of the fingers, that we adjust 
using TGS11 simulations. 
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If the derivatives of the effective gravity are neglected, as well as the horizontal turbu- 
lence, this expression of the thermohaline mixing coefficient is identical to the DIO's one. 
If only horizontal turbulence is int roduced, this expression is similar to that obtained by 



Denissenkov fc Pinsonneaultl (120081 ) 



In the general case, the two limits for the development of thermohaline convection 
are slightly modified. The boundary with dynamical convection, which occurs for large \i- 
gradients, is not exactly R = 1 but Rq = (1 + g e ^)/(l — g e ,r)- The other boundary, for 
small /i-gradients, becomes R = (1 + g e ^)/[r*(l — g e ,r)] instead of R = 1/r . The orders 
of magnitude of these corrections are discussed in the next section. 



4. Orders of magnitude and discussion 

In the previous section, we have introduced the partial derivatives of the effective gravity 
9e,fi — din g e /dha \i and g e ^ = d\ng e /d\nT. We now have to understand their real meaning, 
and derive their orders of magnitude in stellar situations. In the following we first assume 
no external turbulence. We will discuss the effect of adding turbulence later on. 

The effective gravity is g e = g — g R where g R is the radiative acceleration transferred to 
the whole medium by the elements which are individually levitated, namely g R = J2iXig R (i) 
where Xi is the mass fraction of the element i and g R (i) the radiative acceleration acting on 
the same element i. In the following, to evaluate the orders of magnitude of these terms, we 
will simplify the problem by assuming that one major element is the reason for the /i-gradient 
inversion (e.g. iron). Generalisation to several elements will be done numerically in stellar 
evolution codes. 

Even in the case of one levitating element, the situation is not simple as each ion behaves 
in a different way, which is precisely the reason for the selective accumulation. In this case 
we write g R = Xig R {i) where g R {i) represents the contribution to the radiative acceleration 
of all the ions of the considered element, that is: 

X 

g R {%) = EjXijgij = .V,>; ; — ^//,., (22) 

where j represents the various ionisation stages of the element i. 

The variation of g R with /i is a combination of two effects. First the relative mass 
fraction of the heavy element Xi appears as a multiplicating coefficient in g R . Second this 
mass fraction has also an influence on gn(i) through the saturation of the lines. Namely, 
when the abundance of the element increases, the effect of the saturated lines on the total 
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radiative acceleration of the considered element decreases because all the considered ions 
have to share the same available photon flux. 

The variation of /i induced by the accumulation of element i may be written: 

S\Rn^SXi[\-^\ (23) 

where Wj is the average number of particles (including electrons) associated to element % in 
the prevailing physical conditions and Ai its mass number. 

For the considered heavy elements, the second term is small compared to one (otherwise 
there would not be any inverse //-gradient). We neglect it in first approximation and write 

9e,fi as: 

_ <91n g e 1 dg R _ g R (i) Xj dg R (i) . . 

9e '" d In v ~ g dX { g + g dX % ^ 
In this equation, the first term is positive. The second term is negative, difficult to evaluate, 
but its maximum importance is reached if g R {%) varies like 1/Xi , in which case g e>ll vanishes. 
On the contrary, if the saturation effects are neglected, the g e ^ term can reach the order of 
magnitude of g R (i)/g, which may become larger than one (Figure 3). 

The variation of g R with temperature only appears through g R (i) and comes mostly 
from the fact that the various ions of the same elements are not pushed up in the same way 
by the radiative acceleration due to their different atomic structure: 

_ dlng e _Tdg R _T dg R {i) 
9e ' T -dMf- ~g-&T ~ g i ~dT~ (25) 

or: 

g^^lx^xj^-gj^) (26) 



The order of magnitude of the first term may be evaluated with the assumption that g 



varies inside the star approximately like 1/T (IMichaud et al.lll976f ). It is negative, of order 



—g R /g. The second term is more difficult to evaluate, as X^j depends on the Saha equation. 
It may either increase or decrease with T. In the approximation of two ionisation stages, j 
and (j — 1), the fraction of ion j is approximately proportional to exp(—Ej_i/k B T), where 
Ej-% is the ionisation energy of ion (j — 1), so that its contribution to the partial derivative 
of g e may be approximately evaluated as {g R / g){Ej^i/k R T) 1 which may become larger than 
g R /g. In our reference model (1.7 M at 403 Myr), this term never exceeds 15% whereas 
g R /g always remains below 3%. 

It is interesting to note that, due to the signs of these terms, the effect of the radiative 
acceleration on the thermohaline mixing coefficient is often to reduce it, but may also some- 
times increase it. In the regions where the radiative acceleration decreases with increasing 
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radius, the falling blobs suffer less levitation than the surroundings, thereby falling more 
easily. In any case, these terms need to be precisely tested in the computations of stellar 
models. 

Let us now discuss the effect of an external turbulence on the thermohaline convection. 



As already pointed out by iDenissenkov fc Pinsonneaultl ( 120081 ). the presence of horizontal 
turbulence decreases the efficiency of the thermohaline mixing and even suppresses it com- 
pletely as soon as the horizontal turbulent diffusivity is larger than the thermal diffusivity, 
in other words as soon as the horizontal Peclet number is larger than one. In the small 
frame of Figure 1 , the curves represent the thermohaline coefficients computed by including 
horizontal turbulent mixing coefficients of 3000 cm 2 .s~l (dotted curve) and 15000 cm 2 .s~l 
(dotted-dashed curve). These turbulent diffusion coefficients are small, but the consequences 
on the thermohaline convection is important. The basic effect is that the thermohaline co- 
efficient drops for values of Rq decreasing with increasing turbulence. This means that 
thermohaline convection needs a larger /z -gradient to take place. If the external turbulence 
is due to rotational mixing as modelled by lZahnl (119921 ) and coworkers, the importance of the 
horizontal m ixing, much larger than the vertical one, cannot be neglected. As discussed by 
Zahnl (119921 ). whereas the Peclet number is always smaller than one in the vertical direction, 
it may become high and even larger than one in the horizontal direction. In this case, simply 
adding the turbulent diffusion coefficients describing the thermohaline convection and the 
rotation-induced mixing strongly overestimates the global mixing efficiency. 



5. Conclusion 

We have derived a new expression for the thermohaline diffusion coefficient in stars, 
which occurs in case of mean molecular weight inversion, taking into account the influence 
of radiative accelerations and external turbulence. We find that the corrections induced by 
radiative accelerations are small but should not be neglected in the computations of stellar 
models. When relatively abundant elements like iron are strongly pushed up by radiation, 
the momentum transfer to the surroundings leads to a negligible effect on the stellar gas as a 
whole. The only important influence of radiative accelerations on thermohaline convection is 
due to the difference in the radiatively induced momentum inside the fingers (falling blobs) 
and outside (surrounding medium). This difference is small compared to the other effects, 
but in any case the influence of radiative accelerations on the thermohaline efficiency should 
be systematically checked in all computations of thermohaline convection induced by the 
radiative levitation of specific elements. 

These results are important as they confirm the idea that the accumulation of radiatively- 
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levitated elements in stars leads to mixing in a straightforward way. Up to now, atomic 
diffusion including selective radiative accelerations was known to qualitatively account for a 
large number of stellar observations (e.g. chemically peculiar stars). However, atomic dif- 
fusion alone leads to abundance variations much too important for a quantitative fit to the 
abundance determinations. Macroscopic motions like rotation-induced mixing, mass loss, 
etc. had to be invoked to try explaining the observations. Here we show that a fundamental 
process was forgotten in all previous computations: thermohaline instabilities directly in- 
duced by the atomic diffusion process itself. All the computations of radiative levitation of 
chemical elements will have to be revised, taking this effect into account. Note that in the 
case where heavy elements are levitated, helium always sinks due to gravitational settling, 
as the radiative flux is never able to sustain it in its original abundance. This creates a 
stabilizing /z-gradient which opposes the destabilizing one due to the heavy elements. All 
these effects have to be treated together, and the expected result is that heavy elements may 
still accumulate in connexion with the h elium depletion, but not as much as obtained with 



pure diffusion (C.f. iTheado et all (120091 )) . 



Furthermore, these processes have to compete with external turbulence and particularly 
rotation-induced mixing. This will be somewhat complex, as it depends on the ratio of 
the horizontal to the vertical mixing coefficients. We have shown that a large horizontal 
turbulence decreases and may even suppress the efficiency of the thermohaline mixing. On 
the other hand, vertical turbulence increases this efficiency. A combination of all these 
effects, including their close interaction, will have to be studied in real stellar cases. 

SV acknowledges fruitful discussions with Dr. Haili Hu at the Kavli Institute for Theo- 
retical Physics, University of California, Santa Barbara, during the program Asteroseismol- 
ogy in the Space Age, September to December 2011. 
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Fig. 2. — Thermal diffusivity (short-dashed curve) and atomic diffusion coefficients (or par- 
ticle diffusivities) for helium (solid curve) and iron (long-dashed curve) in the 1.7 M , 403 
Myr model as a function of the external fractional mass. 
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Fig. 3. — Radiative acceleration on iron (large frame) and total radiative acceleration on the 
stellar medium (small frame) compared with the local gravity in a 1.7 M Q , 403 Myr model. 
The dashed curve in the large frame represents the gravity, plotted as a function of the 
external fractional mass. The solid curve is the radiative acceleration on iron, which largely 
exceeds the gravity in several stellar layers, leading to the iron accumulation. The dotted 
curve represents the mass fraction of iron in this model, showing the radiatively-induced 
accumulation. The small frame presents, with the same absissa, the ratio of the global 
radiative acceleration (g R = X(Fe)g ra d(Fe)) induced by iron on the stellar medium, to the 
local gravity. It never exceeds 3%. 



